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Abstract. In a companion article (referred hearafter as paper I) a detailed study of the simply 
transitive Spatially Homogeneous (SH) models of class A concerning the existence of a simply transitive 
similarity group has been given. The present work (paper II) continues and completes the above study by 
considering the remaining set of class B models. Following the procedure of paper I we find all SH models 
of class B subjected only to the minimal geometric assumption to admit a proper Homothetic Vector 
Field (HVF). The physical implications of the obtained geometric results are studied by specialising 
our considerations to the case of vacuum and 7— law perfect fluid models. As a result we regain all the 
known exact solutions regarding vacuum and non-tilted perfect fluid models. In the case of tilted fluids 
we find the general self-similar solution for the exceptional type VI_i /g model and we identify it as 
equilibrium point in the corresponding dynamical state space. It is found that this new exact solution 
belongs to the subclass of models n% = 0, is defined for 7 G (§> §) an d although has a five dimensional 
stable manifold there exist always two unstable modes in the restricted state space. Furthermore the 
analysis of the remaining types, guarantees that tilted perfect fluid models of types III, IV, V and VII^ 
cannot admit a proper HVF strongly suggesting that these models either may not be asymptotically 
self-similar (type V) or may be extreme tilted at late times. Finally for each Bianchi type, we give the 
extreme tilted equilibrium points of their state space. 
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1. Introduction 

In a previous article (designated in the following as paper I) we presented a general study of the 
Spatially Homogeneous (SH) models of class A admitting a four dimensional group of homotheties acting 
simply transitively on the spacetime manifold. We concentrated on tilted 7— law perfect fluids and we 
have given various results concerning the existence of an admissible self-similar perfect fluid model in all 
the Bianchi types of class A. As a consequence this analysis has led us to regain the type II solution , 
to show the non-existence of self-similar type VIIo , VIII and IX tilted perfect fluid models 0] and to 
prove the generality of the known perfect fluid solution for type VIo models ^ E] • 

Our intention in the present work is to address the problem regarding the existence of a proper 
Homothetic Vector Field (HVF) for the remaining SH models namely those of class B and to examine the 
physical and dynamical implications of the geometric results with the hope to provide a way of obtaining 
a better understanding of their past and future asymptotic dynamics. 

Section 2 contains some basic elements of the two major frameworks of studying SH cosmologies: 
the orthonormal frame and metric approaches. In particular and adopting the conventions and the 
methodology of 6 , we briefly discuss the structure of the dynamical state space of the SH models, 
coming from the evolution equations and the algebraic restrictions for the case where the spacetime is 
filled with a tilted (in general) 7— law perfect fluid. In section 3, by exploiting the basic relations of the 
metric approach and the results of paper I we find all the (simply) transitively self-similar SH models 

1. e. those models which admit a proper Homothetic Vector Field (HVF) without assuming a specific 
matter content filled the spacetime. From these results we proceed in section 4 and give an analysis of 
their possible physical interpretation. Due to their sound dynamical significance, we apply the general 
geometric results only to the case of vacuum and 7— law perfect fluid models. For illustration purposes 
we reproduce the known vacuum and non-tilted self-similar models of type III, IV, VI;,. In the case of 
tilted perfect fluids we find a new exact solution, representing the general self-similar tilted model of type 
VL^g. Finally section 5 contains a summary of all the known self-similar tilted perfect fluid solutions 
and a brief discussion regarding their importance in the asymptotic behaviour of generic models. 

Throughout the following conventions have been used: spatial frame indices are denoted by lower 
Greek letters a, (3, ... = 1, 2, 3, lower Latin letters denote spacetime indices a, b, ... — 0, 1, 2, 3 and we use 
geometrised units such that 8irG = c = 1. 

2. Preliminaries 

Using the so-called metric approach, one can easily derive the explicit form (in local coordinates) of the 
self-similar metric and the conformally mapped tilted fluid velocity. Due to the fact that the reduced Field 
Equations (FE) have purely algebraic form, the determination of a specific self-similar vacuum or 7— law 
perfect fluid model is straightforward. On the other hand the orthonormal frame approach is used to write 
the FE in an autonomous form and the resulting system of decoupled first order differential equations is 
studied in the dynamical state space of the corresponding model 0. Since every (non-extreme) equilibrium 
point is represented by an exact self-similar solution the mutual use of the approaches is able to show the 
existence or not of a self-similar model and to provide with the necessary tools for the description of their 
future asymptotic dynamics. Therefore and in order the present article to be self-contained as possible 
we find convenient to give the general setup concerning the essence of both approaches which shall be 
used frequently in what follows (the reader is referred to [E| [7| |S1 EH E3 f° r an extensive exposition of the 
formalism used in this section). 

2.1. The metric approach 

In Bianchi models the existence of a Q3 Lie algebra of Killing Vector Fields (KVFs) X Q with 3-dimcnsional 
spacelike orbits C implies the existence of a uniquely defined unit timelike congruence n a (n a n a = — I) 
normal to the spatial foliations C: 

1 9 

n [a;b\ = = n a . b n h <^ -Cn9ab = ™a;b = 0"a& + ~h a b (2.1) 
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where er ao , 6 = n a -bQ 1 , ^ab = 9ab + « Q ^6 are the kinematical quantities associated with the n a according 
to the standard kinematical decomposition of an arbitrary timelike congruence Because n a is 

irrotational and geodesic, there exists a time function t(x a ) such that n a — 8f i.e. each value of t 
essentially represents the hypersurfaces C. Therefore local coordinates {t,x a } can be found in which the 
line element takes the following form: 

ds 2 = -dt 2 + g a0 (t)uj^uj^dx a dx b (2.2) 

where we have employed a basis of invariant vector fields and 1— forms constituting of e Q and its dual 
wf>. 

Furthermore if we assume that the spacetime is filled with a perfect fluid having a linear equation 
of state p = (7 — the Ricci tensor is written: 

-Rah = TM^aMfe H g a b (2.3) 

where fl,p are the energy density and the pressure measured by the observers comoving with the fluid 
velocity u a . The tilted (in general) fluid velocity u a can be decomposed parallel and normal to n a as 
follows: 

u = Tn + TA a uj a = Tn+TB a e a (2.4) 

where B a (t), A a (t) are the frame components of the spatial part of the four- velocity u a and T is a smooth 
function of the time coordinate t satisfying the constraint: 

r = (l-S Q A Q )"5 (2.5) 

and B a A a < 1. 

Finally we note that under the assumption of a TI4 Lie algebra of homotheties, it has been shown in 
paper I that T becomes a constant and both the frame components of the metric and the fluid velocity 
are determined explicitly up to integration constants. 

2.2. The orthonormal frame approach 

Of particular importance in the study of the asymptotic dynamics of SH models, is the reformulation 
of the FE as an autonomous system of first order ordinary differential equations. This has been done 
in where, using the orthonormal frame approach, the FE are written in terms of a set of expansion- 
normalised variables, by using the dimcnsionless time parameter: 

dt 1 dH /H , . 

d-r = H> ^-^ l + ^ H ^ 
where q, H being the deceleration and expansion (Hubble) parameter respectively of the normal timelike 
congruence n a . 

The resulting system consists of the evolution equations for the components E Q/ 3 of the shear tensor 
of the unit normal vector field n a , the spatial curvature A a , n a p of the orbits of the G3 isometry group 
(due to Jacobi identities) and the spatial components of the tilted fluid velocity v a . 

Using the freedom of a time-dependent spatial rotation, we may choose the orthonormal tetrad to 
be the eigenframe of n a /j therefore the contracted form of Jacobi identities n a pAP = implies: 



n a/ 3 






) , A a — A^. (2.7) 



iV 



3 



We note that in types VI?,, and Vllh the evolution equations for n a p and A a can be used to express the 
component A\ in terms of the curvature variables: 

A\ = hN 2 N 3 . (2.8) 
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In addition and following we introduce the shear variables: 



£+ — 2 (^ 22 



£33) : 



2\/3 



(S22 — £33) 



Si 



J_, 

V3' 



1 



J 23, 



1S12 



1 



-13 



(2.9) 
(2.10) 



° V3 ' ° V3 
where for the sake of simplicity we can drop the tilde. 

With these identifications the expansion-normalized variables {£+, Si, £3, £13, N%, A 3 , v a } 
satisfy the set of evolution equations (A.11)-(A.13) and (A. 29) plus the algebraic constraint (A. 18) given 
in [Sj. Consequently the state space is a compact subset V C R 7 (the compactness of V can be shown 
using the generalised Friedmann equation and the inequality Q, > for the energy density parameter). 

We conclude this section by noting that for each Bianchi type, the evolution equations and the 
algebraic constraints can be determined by specialising them in the standard way: 



Type 


N 2 


N 3 


Restrictions 


IV 





A 3 


A 3 > 


V 








None 


VI h 


A 2 


A 3 


A 2 A 3 < 




N 2 


A 3 


A 2 A 3 > 


VI- 1/9 


N 2 


A 3 


A 2 A 3 < 



3. Solution of the symmetry equations 

In this section we will give the general solution of the similarity equations £ng a b = 1^>g a b where 
H = Hn + A Q X Q is the generator of the one-parameter group of homotheties. Using the fact that 
[H,XJ = cH a Xp and [n,X Q ] = it follows H = H(t) and A = A(x a ) where H(t), A(x a ) are smooth 
functions of the spacetime manifold [J. We note that the exact form of the self-similar SH metrics is 
independent from the physical assumption we made for the matter content of the spacetime. However if 
we employ a perfect fluid source for the gravitational field, the well known relation CuRab = £nT a b = 
implies |12j that the fluid velocity is mapped conformally by the HVF i.e. £nu a = ip u a- 

In what follows the coordinate forms of the KVFs, the invariant basis and its dual are those which 
appear in 10 . We recall that in types VI^, and VII^ the group parameters and h arc related according 

to h = — [^t 1 ]-! and h = respectively. 
Type III 

In contrast with the type VI_i/g models (h — —1/9 <^4> Hr = —2 or Kr = —1/2) in which its 
" exceptional" dynamical behaviour do not pass to geometry (at least regarding the determination of the 
corresponding self-similar metrics), type III models (h = — 1 Iir = 0) cannot be treated simultaneously 
within the class of type VI^ metrics due to the "singular" property of the value h,R = in Jacobi identities. 

In local coordinates {t,x, y,z} the KVFs {X Q } and the canonical 1-forms {co a } are: 

Xi = a y , x 2 = a z , x 3 = d x + yd y (3.1) 

lu 1 = e~ x dy, lj 2 = dz, uj 3 = dx. (3.2) 

The non-vanishing structure constants C\ 3 — 1 and the Jacobi identities of the homothetic Lie algebra 
Ha imply that the remaining non-vanishing structure constants C^ 4 are either C| 4 = b or C| 4 = c. We 
consider subcases (in order to pertain the generality of the results we avoid to set some group constants 
equal to 1): 



Case Ai 
HVF 



H = iptd t + D x d x + bzd z 



(3.3) 
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Metric 



A 3 = v 3 t 



Cll t 2< - Dl+ ^/^ Cl2 t( £ 'i+2V'-b)/'A Cl3t (2^+r>i)/V 



5«/3 - Cl2t (D 1+ 2^-b)/4, C22t 2^-b)/4> 



C23 



t (2V-6)/V 



(3.4) 



(3.5) 



Case A 2 
HVF 



Fluid velocity 



Metric 



where 



H = <ptd t + D x d x + D 2 e x d y + cxd z 

Ai = v^^l*, A 2 = v 2 t, 
( 



A, = - 



(3.6) 



+ (3.7) 



5a/3 



521 531 \ 
C22^ 2 .932 



Cnt 2(V+I3i)/V 
C 12 £( 2 V'+£>i)/V> 



c 22 clnt\ c 12 D 2 t^ +D ^l^ 

.932 = t C23 " 



c 2 C22(lnt) 2 2c 2 3clnt 
+ pt-^ ; + c 33 



V> 2 



* 2 . 



(3.8) 



(3.9) 



(3.10) 



Case A 3 



HVF 



H = ijitdt + D 2 e x d v + cxd z 



Fluid velocity 



A, = r, /. A 2 = r_>/. A :i = / ( i; 3 - D2Vl ^ ° V2 lnt 



Metric 



where ,g 33 = t 2 



9a/3 



C33 + 



t 2 (ci 3 



Cllt 2 
Cl 2 t 2 

-D2C11+C12C - 



.921 

c 22 t 2 



531 
532 



lnt) t 2 (c2 3 - g ' c "+ e " c lnf) 533 

P2 c H+ 2 - P 2Cl2C+e 2 C 2 2 (} n f) Z _ 



-,2 _ 2(D 2 c 13 +c 23 c) lnt 



(3.11) 
(3.12) 

(3.13) 
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Type IV 

In type IV models, the KVFs and the dual basis have the following coordinate forms: 

Xi=fy, X 2 = 9 Z , X 9 = d x + (y + z)d v + zd x (3.14) 

lu 1 = er x dy - xe~ x dz, u 2 = e~ x dz, uj 3 = dx. (3.15) 

Since the non-vanishing structure constants are C\ 3 — C 23 — Cf 3 = 1 a similar analysis shows that the 
remaining non-vanishing structure constants C^ 4 of H.4 are C\ A = a, C 24 = a. Therefore we have the 
following results: 



HVF 



H = iptdt + Did x + aydy + azd z 



Fluid velocity 



Metric 



9af3 



( Cll f 2 (£'i+V>-a)/V> 



D\v\ lnt 

521 
522 



(3.16) 

A 3 = v 3 t (3.17) 



where 



y Cl3t (Di+^-a)/^ t (D 1+ 2i,-a)/i, ^ + In ^ ^ 

g 22 = t 2(£>i+</-a)/V> 



2D lCl2 lni £>fcn(lntr 
C22 H ; h 



531 \ 

932 

2 



(3.18) 



(3.19) 



Type V 



In this case the KVFs {X Q } and the canonical 1-forms {w a } are: 
X 1 =d y , X 2 = d z , X 3 = d x + yd y + zd z . 
lu 1 = e~ x dy, uj 2 = e~ x dz, uo 3 = dx 



(3.20) 
(3.21) 



C24 — a. 



It turns out that the structure of the homothetic algebra is C\ 3 = C 23 = 1, C 2 4 = b, 
The analysis of the symmetry equations suggests that we must again consider subcases according to the 
vanishing of the parameter a: 



Case Ax (a ^ 0) 
HVF 



H = iptd t + Did x + bzdy + ayd z 



Fluid velocity 



(D 1 -VlIb + 4,) In t (-0! + ^^+^) In t 

Ai = v\e * + v 2 e * . 
. _ Vi(Ai), t -(£>i + V)Ai 



Metric 



(3.22) 
(3.23) 

v 3 t (3.24) 
a 

2{D 1 -^/ab + 4 > ) lnt 2(r> 1 + v^6 + ^) In t 2(0!+^) 

9ii = ci 2 e * +c 2i e +c 23 i ^ (3.25) 

V , *(fl'll),t - 2 (-Dl + VOSll (0 1 -V^+2^)lnt ( o 1 + vinr+2v)in t 

512 = ■ , 5i3 = ci 3 e +c 3i e (3.26) 
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522 = 



2a 2 



tpt(g 13 ), t - (Di + 2^).gi3 _ +2 

523 — 7 533 — C33t 
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(3.27) 
(3.28) 



Case A 2 



HVF 



H = <ptd t + Did x + bzd y 



Fluid velocity 



Ai=vit * , A 2 = t \v 2 — ) , A 3 = v 3 t 



Metric 



/ Cut 



2CD]-j-jO 



9a/3 = 



(3.29) 



(3.30) 



2CDi+j0 



(ci2-^^) t~ *~ 



C22 

^1 + 2^ 



521 

2fcci2 lnt , b 2 cn(lnt) 2 



4P 



93 
532 



1 \ 



.(3.31) 



Type VI hfl 

The KVFs {X Q } and the canonical 1-forms {cu a } are: 

X 1 =d y , X 2 =d Z} X 3 = d x + yd y + h R zd z . (3.32) 

= e- x dy, uj 2 = e~ hRX dz, uj 3 = dx (3.33) 

where Iir ^= and C^ 3 = C 23 = 1- Hence the remaining non- vanishing structure constants C^ 4 are 
(7^4 = a, C%4 = b. Using these results we get: 



HVF 



Fluid velocity 



Metric 



H = iptd t + Did x + aydy + bzd z 

Cl3t (D 1+ 2^-a)/^ C23t (D lhR +2^-b)/i, ^2 



(3.34) 
(3.35) 

(3.36) 



Type VIl hR 

Finally in type VHh R models we have: 

X 1 =d y , X 2 =d z , X 3 = d x - zd y + (y + h R z)d z . 

uj 1 = [ai - ~Y A ^j d v - A ^ dz -> = A ^ d v + + ~Y A ^j dz > ^ 3 = dx 



(3.37) 
(3.38) 
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where: 



h \ h 

A\ = e "^ x cosira, A2 — e ~ x smwx, w 

w 



(3.39) 



The non-vanishing structure constants are (h 2 R < 4) Cf 3 
expressions: 

C 14 = a, C| 4 = a. 



-C23 = /ir 1 C| 3 = 1 leading to the following 



(3.40) 



HVF 



H = iptdt + Did x + aydy + azd z 



Fluid velocity 



V\2 COS 



Dili), \ . f D lW ^ \ 

— ■ — in t + ii 2 i sin — ■ — In t 
tp J \ ip J 



ipt (Ai) jt + (o-^)A x 



A 3 = u 3 t. 



Metric 



5ii 



512 



4*/* {wDiE^pj cos (^f^) + 4S 2W 2 pi sin (^^) } 



64D^ 5 pf 



^ (5ll) f + 2(0- VOfl 1 !! 



.913 = £ 



2L>i 

(pi+2V)/2V 



DiwlnA / Diwlnt 

ci 3 cos I ■ + c 3 i sin 



.922 



523 



V^ 2 (ffn) « - ^ (511) t (Pi + ^ - 2a) + 2 £>? + (tp - a) Pl \ 511 



tpt (gi 3 ) t + (a- 2tp)g 



13 



2Z^ 



533 = C 3 3* 



(3.41) 

(3.42) 
(3.43) 

(3.44) 
(3.45) 
(3.46) 

(3.47) 
(3.48) 



where E±, E2 are constants and we have set p\ = DiIir + 2 (tp — a). 
4. Exact Solutions 

The well known feature of self-similarity to reduce the FE to a purely algebraic form, makes the 
determination of vacuum and perfect fluid models straightforward. Following the methodology of paper I 
and for illustration purposes, we reproduce all the known vacuum and non-tilted perfect fluid solutions. In 
addition we find the general tilted perfect fluid model for the exceptional class VI_i/g. For completeness 
we also give, for each Bianchi type, the extreme! tilted equilibrium points of their state space. 

Vacuum plane wave models of type III 

This exact vacuum solution can be found from the case Ai of section 2 and has been given in |13|. 
The HVF is H = iptdt + D\d x + bzd z and the constants c a p, tp, b are: 



C13 = C23 = 0, en = 



c\ 2 {Di-b) 



C33 



2c 22 £>l ~ (£>! + ft)" 

\ We recall that an extreme tilted equilibrium point satisfies v a v a = 1. 



tp = b 



(4.1) 



Self-similar Bianchi models: II. Class B models 



9 



where without loss of generality we can set 6=1. 

We note that the above spacetime is algebraically special (Petrov type N) since it admits the 
gradient null KVF l a — e -&>yC D i+ 6 ) (pt + Dl b ^ b d x ) which essentially represents the repeated principal 
null direction of the Weyl tensor. 

Non-tilted perfect fluid models of type III 

It is a special case of the non-tilted perfect fluid solution found by Collins jTl] for h = — 1. The HVF 
is H = xptdt + bzd z and the constants c a 0, ip, 6 are given by: 



Cl3 = C 2 3 = C12 = 0, c\x = c 22 = 1 



(4.2) 



C33 



(2- 7 ) (3 7 -2)' 



4(1-7) , trf 



Tilted perfect fluid models of type III 

Using the results of section 2 we can show that there is no physically acceptable tilted perfect solution 
of type III. We note that this conclusion can be also confirmed from the set of the tilted equilibrium points 
of the dynamical state space given in |15| . 



Vacuum plane wave models of type IV 

From the self-similar metrics lf3*.l£|) we can find the general type IV vacuum solution, first given in 
|lb|. In this case the HVF is H = tptdt + D\d x + ayd y + azd z or, by performing a suitable change of the 
basis of the homothetic Lie algebra, H = tptdt + (D± — a) d x — azd y . In addition the constants c a p, ip are 
given by: 



Di( c? 1 -4c? 2 ) + 4ac? 2 

4c 11 (a-D 1 ) ' " M D\ 



ci 3 = c 23 = 0, c 22 = ^ ;; 12 ^ 12 , 033 = ^, tf = a. (4.3) 



It can be verified that the above one parameter family of models is also algebraically special and admits 
the gradient null KVF l a = e~ ax ^ Dl (d t + %d x ). 

Non-tilted and tilted perfect fluid models of type IV 

It is found that no orthogonal or tilted perfect fluid self-similar models of type IV exist. Nevertheless 
there exists the following extreme tilted equilibrium point 17 : 

N 2 = 0, N 3 = N 3 , v a v a = 1, 

V ^ 1 A 

2ji = , vi — I, Ai = 

6 ' 2 

„ 3q(2-q)-Ni 
6 



Models of type V 

The analysis shown that no self-similar vacuum or perfect fluid models exist for type V. This result 
can also be proved by using the set of evolution equations given in [5] (see also JS]). However there exists 
the following extreme tilted equilibrium point ^1 ^] : 

•^2 — N$ — 0, v a v a = 1, 
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E 2 = (Ai-l) 2 , S3-Si 3 -S 1 = S_=0, £+ = -§, 

Wi = l, A x = ^— — , n = 2A 1 (l-A 1 ), 0<7<2. 



Vacuum plane wave models of type VI/j R 

In this case the HVF is H = iptdt + D\d x + ayd y + bzd z and the constants c a /3, ip are given by |13| : 

c 2 2 + 1) [Di + 1) - a - b] 



Cl3 = C 2 3 = 0, Cn 



2c 22 [-Di + 1) - a - bh R ] 

(4.4) 



(ah R - b) 2 a/ifl. - 6 

[£>i (h R - 1) + a - 6] 2 hK-1 



The gradient null KVF is 



"r>l(h R -l) + a-b 



9 t + 



Oi(/in-l) + a-b 



9 3 



y * T (ah R -b)t 

As we have mentioned in section 2, in order to determine the corresponding self-similar solution for 
the exceptional type VI_i/g models we can use the geometric results (|3. 3411 — (f 3 . 3611 and solve the FE by 
setting h R = —2 or h R — —1/2 (for these values of the group parameter the 01 or 02 component of 12.311 
vanish and essentially represents the "exceptional" behaviour of the VI_i/g models in terms of the FE). 
The resulting solution has been first given by Robinson and Trautman (RT) |2U] (see also ^3]) an d has 
the following form (C13 ^ 0): 



n 75 

C12 = C23 = 0, C33 = — , 



24c 2 3 
125 : 



Di = — , 
3' 



v 6 ' 



(4.5) 



Non-tilted perfect fluid models of type ~Vlh R 

The family of non-tilted perfect fluid models has been found by Collins |14| . The HVF is given by 
equation H3.34JI and the constants c a p,i/j,a,b are: 

n , 7 2 (hR I) 2 

C13 = C23 = C12 = 0, cn = c 22 = 1, C33 = 



(2- 7 ) (3 7 -2) 

4[^(7-l) + fe fl7 + 7-l] , _ j(Di-a) (h R - 1) 
M ~ 7 2 (^-l) 2 t 2 ' M2-37) ( ] 

Di (h% + 1) - a 
h R 

The exceptional type VLx/9 non-tilted perfect fluid solution can be found similarily by setting h R = — 2 
(or h R = —1/2) and has been given in |21| : 

75 56 

C12 = c 23 = 0, c 33 = — , ^ = yj Di = 0, a = 2b 

(4.7) 

_ 9(125cn - 24^) _ 10 
^ ~ 25(75cn - 8c 2 3 )i 2 ' 7_ T' 
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Tilted perfect fluid models of type VI/ lR 

It is convenient to employ the constants pi,P2- 

Di=a+pii)-2il), a = bh R ~iP{ Pl +p 2 -A). (4.8) 
The frame components of the metric and the fluid velocity become: 

9a/3 = c 12 t^ h R~ lS > +hB ^ 2 ~ p ^ +p1 ^/^ c 2 2i 2[h(/l «" 1)+ ' ifl ' / ' (2 " P2)+ ' / ' 1/ ' / ' 532 (4-9) 
V " Ci 3 i Pl C23t [b(hl~l) + h R ,p(2- P2 )+24,]/,p C33t 2 J 

Ai = wii" 1-1 ' A 2 = V2i \Kh\-i)+h R W- P2 )+Wi> ^ As = V3t ( 4 1Q ) 

whereas the HVF assumes the form: 

H = i>td t + [bh R -i>(p 2 - 2)] d x + [bh R - V>(pi +p 2 - 4)] yd y + bzd z . (4.11) 

For the case h R ^ —2 (or h R ^ —1/2) the general self-similar solution has been found recently |15|. Here 
we give the corresponding general solution for the type VI_i/ g models. 
We define the parameter s according to: 

7 = 2 . (4.12) 
' 2s + 1 v ' 

Then, the 00-conservation equation implies that the constant p\ is given by: 

ft = 2(fo+ 4 "- K . (4.13) 
The various integration constants are: 

ci 2 = ci 3 = V! = 0, V = t; 4.14 

P2 - 2 



C22 



C33 



c 2 3 (p 2 -2) 2 [2(3s-l)-p 2 ][-9p 2 +4p 2 (3-2s)+624s 2 -208s+12] 
96s[25p|-20p 2 (6s + l)+4(6s+l) 2 J 

48[3p 2 +4p 2 (8s-3)-4(12s 2 +8s-3)] 
(p 2 -2) 2 [9p 2 +4p 2 (2s-3)-4(156s 2 -52s+3] 

s(5p 2 -12s-2)[2(3s-l)-p 2 ][-9p 2 +4p 2 (3-2s)+624s 2 -208s+12] (4-15) 
2[p 2 +2(2s-l)][p 2 (31s-6)+4p 2 s(78s-17)-4(1476s 3 -660s 2 + 101s-6)J 

36(p 2 +4s-2)[p 2 +2(2s-l)] 



9p|+2p 2 (4s-15)-12p 2 (52s 2 -16s-l)+8(156s 2 -52s+3) 



_ 3c 23 [p 2 +4p 2 (s-l)-4(2s-l)] 
V2 ~ 4[2(6s+l)-5p 2 ] ' 



The constant p 2 is related with the "state parameter" s according to: 

2[-6V36s 4 -204s :i + 133s 2 -22s+l|4s-l|+s(336s 2 -214s+19)] 
P 2 ~ 160s 2 -91s+6 



and the energy density of the model is: 



_ _ (2s+l)[p 2 (31s-6)+4p 2 s(78s-17)-4(1476s 3 -660s 2 + 101s-6] 
f 1 — 144s(p 2 -12s+2)t 2 



(4-16) 



(4.17) 



In order to ensure the positivity of the energy density and the real values of the constant P 2 (together 
with the overall signature of the metric) this new exact solution is only defined when: 



6'4/ ' V 3' 2 
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It is interesting to note that, similarly with the types VIo and VI^, the above rotating cosmological model 
admits the hypersurface orthogonal KVF X 2 = d z and belongs to the subclass n" = 0. In fact the present 
solution arises as an equilibrium point of the tilted perfect fluid Bianchi type VT_x/g dynamical state 
space which we denote it as £ and has the following kincmatical and dynamical quantities: 



Ar _ / 3[g(37-6)+77-6][g 2 (47 2 -47-12)-2q(37 2 -227+12)-377 2 +487-12] 

iV 3 - - V 167 2 l 9 (7-2)+2(7-l)J 

N 2 = -N 3 , £_ =0, £+ = -|, v 2 = -v 3 



v^3[g(57-12)+2(77-6)] ^ _ V6(2-7)[2 ? 2 ( 7 -3)+3 g (37-4)+77-6] 
67 ' 3 12 7 V(9+l)[?(7-2)+2(7-l)] 



„ _ V6(2-7)[2 ? 3 (7-3)+ ? 2 (ll7-18)+2g(87-9)+77-6] 
13 12 7 V(«+l) 3 [?(7-2)+2(7-l)] 



n 



3(2 g -3 7 +2) _ 3^/2(2- 7 )(g+l)(2g-37+2) 2 9(2q-3 7 +2) 2 [3g( 7 -2)+77-6] 

2^37 ' 3 87V 3 7V9(7-2)+ 2 (7- 1 ) ' l&N^ [q( 7 -2)+2( 7 -l)] 

_ [9(37-6) + 77-6][2g( 7 -3)+77-6] [4q 2 (4- 7 )+2<;(8- 15 7 )+277 2 -267] 



127 2 (2 9 -3 7 +2)[g(7-2)+2( 7 -l)] 



V 2 _ [3<?(7-2)+77-6][4g 2 (27 2 -87+9)+4g(97 2 -297+18)+497 2 -847+36] 
^ ~~ 12 7 2 [g(7-2)+2(7-l)] 



_ 7-v/(7- 1 )( 73 7 3 -2537 2 +2407-36)|37-4|-2(7-l)(337 3 -1797 2 +2947-144) 
1 ~ 6(7-1) (7-2) (37 2 - 197+24) 

We also note that there exists the following extreme tilted equilibrium point |19j : 



N 2 = N 2 , N 3 = N 3 , v a v a = 1, 



^ 2 _ 3(7V 3 - N 2 ) 2 - 4(N 2 N 3 + 6V-iV 2 iV3 - 9) „ _ v „ _ n 

^ — ^7 > ^3 — ^13 — ^- — U, 

36 

„ ^(iVa-^) „ 3 2(3 - y/-N 2 N 3 ) 

Si = ~ 6 , S + = -2' « = § , « 1 = 1, 

3(iV 3 - 7V 2 ) 2 - 4(N 2 N 3 + 37^1^) _ _ 
" = 7^ j < 7 < 2. 



Models of type VII, lJ? 

The question of the existence of tilted perfect type Vllh models is notoriously difficult to answer due 
to the complexity of the self-similar metrics (I3.44|) - (I3.48[) . In fact the resulting FE are difficult to handle 
analytically even to the case of expressing the known Vllh vacuum solution (see e.g. [7] page 192 for 
an elegant form of this solution and [22] for a complete study of its stability properties against vorticity, 
shear and Weyl curvature perturbations). Nevertheless, using the set of evolution equations we can show 
that the assumption of the existence of a proper HVF is incompatible with tilted perfect fluid models 
and that type Vllh models possess the following extreme tilted equilibrium point (this equilibrium point 
has been also given in |17|): 

N 2 = N 2 , N 3 = N 3 , v a v a = 1, 



vj (JV 3 - N 2 ) 2 + \2{hN 2 N 3 - 2VhN 2 N 3 + 1) 



s2 = v_s fz_ v L j S 3 = E 13 = S _ = 0, 
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Si 




|, q = 2(1 - y/hN 2 N 3 ), Vl = l, 



n = 



(^3 - ^2) 2 + 12(/tA 2 A 3 - VhN 2 N 3 ) 
6 



< 7 < 2. 



5. Discussion 

The main concern of papers I and II can be regarded having two main branches: 

a) a geometric nature concerning the determination of all the SH geometries restricted only by the 
requirement of admitting a proper HVF and, 

b) the physical implications of the general geometric results by finding, whenever they exist, the 
corresponding vacuum or perfect fluid models. 

The underlying importance of a) is that we have not incorporated a specific form for the matter 
fields filled the spacetime therefore we can associate the self-similar metrics with more general matter 
configurations with a view to analyse their physical significance. On the other hand the well established 
aspect of the SH equilibrium points as past or future attractors for general vacuum or perfect fluid models 
necessitates the knowledge of all the self-similar equilibrium points in order to gain deeper insight into 
their asymptotic dynamics. This was accomplished in the series of articles I, II, [31 0] m which the 
complete set of tilted perfect fluid self-similar models have been found. For convenience we collect them 
in table 1 together with the corresponding references. 

As an immediate application of the geometric results of section 3, we reproduced the known exact 
vacuum and non-tilted solutions for the Bianchi types of class B and we have found the general form 
of the self-similar tilted type VI_i/g models. As we have seen in section 4 this new exact solution arise 
as equilibrium point in the VI_i/g state space and shares many of the common properties of type VIo 
and VI/j (h ^ —1/9) models and important differences as well. For example the stability analysis of 
vacuum and non-tilted equilibrium points |25j has shown that the Collins VI_i/g solution is stable for 
7 € (§, ip). At the arc of equilibrium points 7 = ^ 24 there is an exchange of stability with the RT 
vacuum solution jSJEOj which is stable whenever 76 (-^, |). As a consequence we expected that tilted 
models will be future dominated against the RT model when 7 > |. Indeed the equilibrium point E is 
found to have a five dimensional stable manifold for 7 € (§>§)• However it can be shown that, at least 
for the case where = 0, there exist always two unstable modes in the restricted state space. One can 

made a step further and conjecture that for the type VI/, and VL1/9 models the values 7 = 2 ^^? - 

and 7 = | respectively represent line bifurcations similarly to the type VIo models with 7 = | . However 
no such bifurcations are found in both models. 

Furthermore, for the remaining types, the mutual use of the metric and orthonormal frame 
approaches guarantees that tilted perfect fluid models of types III, IV, V and VII^ cannot admit a 
proper HVF indicating that these models either may not be asymptotically self-similar (a preliminary 
analysis suggests that this behaviour occurs in type V models) or may be extreme tilted at late times 
(see |17) for a thorough discussion on the asymptotic dynamics of type IV models). 

We conclude by noting that an interesting aspect of the analysis presented in papers I and II is 
that vacuum and tilted perfect fluid SH self-similar models exhibit further geometric constraints coming 
from the existence of covariantly constant null vector fields or a hypersurface orthogonal KVF. Therefore 
we believe that it will be interesting to study how the assumptions of the self-similarity and a specific 
dynamic description lead to extra restrictions on the geometry of the SH models. 
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Table 1. This table contains all the self-similar tilted perfect fluid SH models and the relevant range 
of the state parameter 7. 



Type Self-similar Model Reference 



i $ EH1 

II 7 G (f, 2) 2 

VIo 7 6[f,f) HE] 

viio $ n 

VIII $ 

IX 

in $ US] 

IV $ paper II 

V $ paper II 



vih U51 



VU h $ paper II 

VI-i/9 7 €(§,§) paper II 
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